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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  MIIRIC  (SI) 
UNITS  OF  MLASMRkMhNT 


U.S.  customary  units  of  measurement  used  in  this  report  can  he  converted 
to  metric  (SI)  units  as  follows: 


l ri  — ■■  <n  -•» 

Multiply 

by 

To  obtain 

inches 

25.4  " 

millimeters 

2.54 

cent imeters 

square  inches 

6.452 

square  cent imeters 

cubic  inches 

16.  39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  motors 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 10" 

3 kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

ki lograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

^i'o  obtain  Celsius 

(C)  temperature 

readings  from  Fahrenheit  (F)  readings. 

use  formula:  C = 

(5/9)  (F  -32). 

To  obtain  Kelvin 

(K)  readings,  use 

formula:  K = (5/9)  (F  -32)  + 273.15. 
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SYMBOLS  AND  DEFINITIONS 


t t irne 


ox  horizontal  axis  at  S K I.  parallel  to  the  (initial)  beach  profile 


oy  horizontal  axis  at  S W I.  perpendicular  to  the  (initial)  beach 
prof i le 


l)  beach  depth  (depth  beyond  which  sediment  transport  is  negligible) 


a wave  angle  with  beach  profile 


aQ  wave  angle  with  beach  profile  at  infinity 


Q longshore  transport  (littoral  drift)  discharge 


constant  = tt 


i 

D da 


u parameter  u 


E (u)  Fresnel  integral  = E(u) 


A.  /V 

vu  Jo 


l length  of  groin 


tj  time  for  the  beach  profile  to  reach  the  end  of  the  groin 


t'  transform  time  t ' = 0.62tj 


B sinusoidal  beach  amplitude  (at  time  t = o) 


■ 


parameter  related  to  beach  wavelength  1.:*  = I ' 1 K 


parametric  value  of  x defining  volume  of  beach  dumping 

parametric  value  of  y defining  volume  of  beach  dumping 

parameter  used  to  define  hypocycloid  beach  profile  between 
headlands 


coefficients  used  in  the  littoral  drift  formula  to 
characterize  the  effect  of  wave  angle 


breaking  wave  height 


water  depth  at  inception  of  wave  breaking 


group  velocity 


littoral  drift  constant  6.42  x 10  ' 

distance  of  shoreline  from  a horizontal  axis  parallel  to  the 
initial  beach  profile 


distance  of  the  offshore  beach  limit  from  a horizontal  axis 
parallel  to  the  initial  beach  profile 


■> 


v 


1 


equilibrium  distance  y 


onshore-offshore  transport  per  unit  length  of  beach 
onshore-offshore  transport  parameter  (dimension  1.1  *) 
longshore  sand  transport  discharge  in  shallow  water 
longshore  sand  transport  in  deeper  water 

V1  - y2 

distance  of  the  beach  profile  to  a spiral  center 
angle  parameter  in  mathematical  description  of  hooked  bays 
spiral  angle  in  mathematical  description  of  hooked  bays 
depth  of  hooked  bays 


distance  between  headlands 


_ z /vor 

r/-£cediy^  vbpa 

MATHEMATICAL  MODI:  I . I NCI  01  SHORELINE  EVOLUTION 

by 

Bernard  Le  Mehaute  and  Mills  Soldate 

1.  INTRODUCTION 

This  interim  report  presents  a critical  literature  survey  on  the 
subject  of  mathematical  modeling  of  shoreline  evolution.  Hopefully, 
this  review  will  lead  the  way  in  establishing  a flexible  and  practical 
numerical  method  suitable  to  predict  shoreline  evolution,  resulting 
from  the  construction  of  navigation  and  shore  protection  structures  in 
the  Great  Lakes. 

To  focus  attention  on  the  most  pertinent  literature,  the  subject 
under  consideration  is  limited  to  long-term  shoreline  evolution  as 
defined  below. 

Three  time  scales  of  shoreline  evolution  can  be  distinguished: 

(a)  Geological  evolution  taking  place  over  centuries; 

(b)  long-term  evolution  from  year-to-year  or  decade;  and 

(c)  short-term  or  seasonal  evolution  and  evolution 
taking  place  during  a major  storm. 

Associated  with  these  time  scales  arc  distances  or  ranges  of  influ- 
ence over  which  changes  occur.  The  geological  time  scale  deals,  for 
instance,  with  the  entire  area  of  the  Great  Lakes.  The  long-term 
evolution  deals  with  a more  limited  stretch  of  shoreline  and  range  of 
influence;  e.g.,  between  two  headlands  or  between  two  harbor  entrances. 

The  short-term  evolution  deals  with  the  intricacies  of  the  surf  cone 
circulation;  e.g.,  summer  profile-winter  profile,  bar,  rhythmic  beach 
patterns,  etc. 

Lor  the  problem  under  consideration,  long-term  evolution  is  of  pri- 
mary importance,  the  short-term  evolution  appearing  as  a superimposed  per- 
turbation on  the  general  beach  profile.  Evolution  of  the  coastline  is 
characterized  by  low  monotone  variations  or  trends  on  which  are  super- 
imposed short  bursts  of  rapid  development  associated  with  storms. 

The  primary  cause  of  long-term  evolution  is  water  waves  or  wave- 
generated  currents.  Three  phenomena  intervene  in  the  action  which 
waves  have  on  shoreline  evolution: 

(a)  Erosion  of  beach  material  by  short  period  seas  versus 
accretion  by  longer  period  swells; 
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lb)  effect  of  (lake)  level  changes  on  erosion;  and 
(c)  effect  of  breakwaters,  groins,  and  other  structures. 


Even  though  mathematical  modeling  of  shoreline  evolution  has  in- 
spired some  research,  it  has  received  only  limited  attention  from 
practicing  engineers.  The  present  methodology  is  based  mainly  on 

(a)  the  local  experience  of  engineers  who  have  a deep  knowledge 
of  their  sectors,  understand  littoral  process,  and  have  an  inherent 
intuition  of  what  should  happen;  and 

(b)  movable-bed  scale  models  that  require  extensive  field  data 
for  their  calibration. 

In  the  past,  theorists  have  been  dealing  with  idealized  situations, 
rarely  encountered  in  engineering  practice.  It  seems  that  mathematical 
modelers  have  long  been  discouraged  by  the  inherent  complexity  of  the 
phenomena  encountered  in  coastal  morphology.  The  lack  of  well-accepted 
laws  of  sediment  transport,  offshore-onshore  movement,  and  poor  wave 
climate  statistics  have  made  the  task  of  calibrating  mathematical 
models  very  difficult. 

Considering,  on  one  hand,  the  importance  of  the  subject  of  deter- 
mining  the  effect  of  construction  of  long  groins  and  navigation 
structures  and  on  the  other,  the  progress  which  lias  been  made  in 
determining  wave  climate  and  littoral  drift,  it  now  appears  that  a 
| mathematical  approach  could  be  useful. 

I The  complexity  of  beach  phenomena  could,  to  a large  extent,  be 

taken  into  account  by  means  of  numerical  mathematical  scheme,  (instead 
of  in  closed-form  solutions),  dividing  space  and  time  intervals  into 
small  elements,  in  which  the  inherent  complexity  of  the  morphology 
could  be  taken  into  account. 

| Furthermore,  better  knowledge  of  the  wave  climate,  a necessary  in- 

put, will  allow  a better  calibration  of  coastal  constants  such  as 
found  in  the  littoral  drift  formula. 

This  study  emphasizes  the  relative  importance  of  various  reports  and 
reviews  the  most  important  ones.  Conclusions  based  on  this  review  arc 
presented,  pointing  out  the  deficiencies  of  the  state-of-the-art.  (Sub- 
i sequent  investigators  should  attempt  to  bridge  the  remaining  gaps.) 

The  reports  arc  presented  individually,  primarily  in  chronological 
order.  Two  milestone  developments  from  this  survey  are  reports  by 
Pel nard-Considere  (1956)  and  by  Bakkcr  (1968b).  Others  arc  extensions 
and  refinements,  experimental  verifications,  support  papers,  numerical 
procedures,  and  side  issues,  including  the  latest  developments  on 
"hooked  beaches"  or  crenulate-shaped  bays. 
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11.  nil:  FIRST  MODEL  ( PELNARD-CONSIDERE ) 


The  idea  of  mathematically  formulating  shoreline  evolution  is  attri- 
buted by  Bakker  (1968a)  to  Bossen,  but  no  reference  to  Bossen  is  given. 
The  first  report  which  appears  in  the  literature,  on  mathematical  model- 
ing of  shoreline  evolution,  is  by  Pelnard-Cons iderc  (1956).  His 
theoretical  developments  were  substantiated  by  laboratory  experiments 
made  at  Sogreah  (Grenoble),  France.  The  experimental  results  fit  the 
theoretical  results  very  well.  It  is  surprising  that  such  relatively 
s i mp 1 e theory  has  not  been  more  frequently  applied  to  prototype  cases  by 
the  profession  (at  least  as  it  would  appear  from  the  open  literature),  a 
fact  which  may  be  attributed  to  the  lack  of  knowledge  of  wave  climates. 


Pelnard-Considere  assumed  that: 

(a)  The  beach  profile  remains  similar  and  determined  by 
the  equilibrium  profile.  Therefore,  all  contour  lines  arc 
parallel.  This  assumption  permits  him  to  consider  the  problem  to 
be  solved  for  one  contour  line  only. 

(b)  The  wave  direction  is  constant  and  makes  a small  angle 
with  the  shoreline  (<20°). 


(c)  The  longshore  transport,  Q , is  linearly  related  to  the 
tangent  of  the  angle  of  incidence  a • (Q  = f(a) , f(a)  = tan  a). 

(d)  The  beach  has  a fixed  ( i 1 1 -def  i ned ) depth,  D (F;ig.  1). 

II  is  a factor  relating  erosion  retreat  to  volume  removed  from 
profile,  which  could  be  defined  by  the  threshold  velocity  of 
sand  under  wave  action.  A practical  method  of  determination 

of  I)  is  given  in  Section  V 1 1 1 . 

Despite  the  crudeness  of  these  approximations,  the  Pelnard-Considere 
model  can  be  considered  as  a milestone  in  demonstrating  the  feasibility 
of  mathematical  modeling  of  long-term  shoreline  evolution.  For  this 
reason,  it  is  judged  useful  to  describe  in  some  detail  his  theoretical 
development . 

Consider  an  axis,  ox  , parallel  to  the  main  coastal  direction  and 
an  axis,  oy  , perpendicular  seawards  (Fig.  2).  The  angle  the  deepwater 
wave  makes  with  the  axis,  ox  , is  a . The  angle  of  the  wave  with  the 
shoreline  a at  any  location  is  assumed  to  be  small;  therefore, 


- 1 3y 

a = a - tan 

o 3 x 


3y 

a - r^-  or 
o 3x 


3 v 

a - a = - —~ 
o 3x 


(1) 


(y  = f(x,t)  gives  the  form  of  the  shoreline  as  function  of  time  t). 

The  littoral  drift  Q is  a function  of  angle  incidence  a and  can  be 
put  into  a Taylor  series: 
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(2) 


u = ct 

o 

in  which  Q denotes  the  transport,  Q , when  the  angle  of  the  wave 
incidence  is  aQ  . Substituting  equation  (1J  into  equation  (2)  yields: 

«■<*,-  [s 

During  the  interval  of  time,  dt  , the  shoreline  recedes  (or  accretes)  by 
a quantity  dy  . Therefore,  the  volume  of  sand  which  is  removed  (or 
deposited)  over  a length  of  beach,  dx  , is  I)  dx  dy  . The  quantity  is 
equal  to  the  difference  of  longshore  transport  during  time,  dt  , between 
x and  x + dx ; i . e . , 

Q dt  and  : and  ( Q + dx  ) dt  ; 

d X 


1 

a = a I 3x 

o J 


(3) 


i.e.  , 


57  -* 


Therefore , 


U dx  dy  = 


dxdt 

3x 


or 


3y  _ l 3Q 
3t  D 8x 


(4) 


Substituting  the  expression  for  Q , o being  small,  and  defining 


K _ I dQ 
D da 


yield: 


K 


(6) 


which  is  the  well-known  diffusion  or  heat-flow  equation. 

K is  approximately  constant  at  a given  site.  Bakker  (1968a)  found 

K equal  to  0.4  x 106  cubic  meters  per  meter  depth  per  year,  at  an  exposed 
site  along  the  coast  of  the  Netherlands.  Equation  (6)  demonstrates  that  the 

rate  of  accretion  or  (erosion),  , is  linearly  related  to  the  curvature  of 

dt 


15 


the  coast,  the  derivative  of  the  longshore 
to  the  angle  of  the  wave  incidence,  I 


a 


tional  to  the  beach  depth,  D . 


transport  rate  with  respect 
, and  inversely  propor- 

= a 

o 


The  above  equation  will  be  recognized  as  the  well-known  diffusion 
equation.  A number  of  classical  solutions  of  mathematical  physics  are 
applicable  to  the  diffusion  equation  when  boundary  conditions  are 
specified.  Pelnard-Considere  (1956J  applied  his  theory  to  the  case  of 
a littoral  barrier  or  long  groin.  This  case  is  reviewed  below; 

The  longshore  transport  rate  along  a straight,  long  beach  is  sudden 
ly  stopped  by  the  construction  of  a long  groin  built  p'crpcndicular  to 
the  beach  (see  Fig.  2).  The  boundary  conditions  are: 


(a)  y = o for  all  x when  t = o which  characterizes  an 
initial  straight  shoreline. 

(b)  At  the  groin,  the  longshore  transport  rate  Q = o which  is 
realized  when  the  waves  approach  the  shore  normally;  i.c.,  when 


( 


i 


3y  = -tan  a at  x = o , 

— — o 

3x 


(c)  9y_  = o at  a large  distance  updrift  ( x -*  °° 

9x 


and  Q = Q 

a 


Q is  the  steady-state  longshore  transport  along  a straight  beach 
for  tRe  given  wave  conditions.  The  solution  for  the  given  boundary 
conditions  is: 


tan 

a 

* 

- 

y - — 

o 

V4Kt  exp  (u“)  - x V*"1  E 

(u) 

' (73 

\A  | 

- 

- 

where  u = ^ , 

and 

E (u)  is  the  Fresnel  integral 

> 

E (u)  = 

2 

f 2 

1 e~u  du 

■\Ztt 

■'ll 

(8) 

Values  of  E (u) 

or  more  frequently,  4>  (u)  = 1 - E 

(u) 

, can  be  found 

in  tabulated  form  as  given  in  Table  1: 
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Table  1.  u versus  0 (u) . 


u 0.1  0.2  0.3  0.4  0.5  0.7  0.9  1.2 


0 (u)  0.112  0.223  0.328  0.428  0.520  0.667  0.796  0.910  0.995 


Fig.  2 illustrates  the  shoreline  evolution  as  defined  by  equation  (8). 
It  is  interesting  that  these  curves  are  horaothetic  with  respect  to  the 
origin  o ; i . e . , 


The  horizontal  lengths  grow  with  t , and  in  particular, 


oy  = 


2 VkT 


A tangent  to  the  shoreline  at  the  groin  intersects  the  initial  shoreline 

defined  by  y = o at  a point  a distance  of  2 VlCt/m  updri ft  from 
the  groin. 

The  ratio  of  the  area  of  sand  accumulation,  such  as  is  in  oyxD  , to 

D 

the  area  of  sand  contained  in  the  triangular  fillet,  oyx  , is  1.56  and 
the  distance  ox^  a 2.7  ox  . This  ratio  permits  rapid  assessment  of 

the  total  amount  of  sand  accumulated  updrift  from  a single  measurement 

of  the  angle  a , and  determination  of  D as  shown  in  Section  IV. 
b o 

The  end  of  the  groin  of  length,  oy  = l , is  reached  when 
2 

t = t = 1 71  • (10) 

* 4K  tan  a 


When  t t j , the  boundary  conditions  must  be  modified  since  the  groin 
no  longer  traps  all  the  sand  but  bypasses  some  of  it. 

If  the  same  theory  is  applied  to  the  beach  downdrift  of  the  groin 
and  if  assumed  that  the  wave  diffraction  effects  are  negligible,  the 
beach  is  eroded  in  a form  symmetric  with  the  updrift  accretion. 


: 

r 


When  t = t , the  end  of  the  groin  is  reached  hy  the  shoreline  and 
sand  begins  to  ^ be  bypassed  around  the  groin. 

The  boundary  condition  at  the  groin  becomes  oy  = St  (constant)  for 
t > t . The  solution  then  becomes  (Fig.  3): 


The  curves  representing  the  shoreline  become  homothetic  with  respect  to 
the  axis  oy  ; i.e., 


AB  _ AC 
A ' B ' AT' 


etc . 


The  area  between  the  shoreline  and  the  ox  axis  (°yox  ) is  8ivcn  ’’>  • 


The  area  of  triangular  fillet,  °>'0X0  > 


St 


N/nKt 


Hence , 


1 = i = 1.27  (12) 

TT 

I 

and 

ox'  = 2x 

o 

The  shoreline  as  described  by  equation  (7)  at  time  t = tj  is  slightly 
different  from  the  shoreline  defined  by  equation  (11)  at  t = tj  as 
shown  in  Fig.  4. 


oy  x 
1 o 

oy  x 
1 o o 


The  volume  of  sand  defined  by  both  curves  is  equal  when  the  time  tj  of 

equation  (7)  is  replaced  by  the  time  t'  in  equation  (11)  in  such  a way 
that 


71 

To 


= -nr  i-e-> 


= 0 . 62t 


(13) 
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Figure  3.  Successive  beach  profiles  updrift  of  a groin  after  sand 
bypassing  (from  Le  Mehaute  and  Brebner , 1961). 


p 


Therefore,  the  shoreline  evolves  initially  as  represented  by  equation 
(7).  Then,  when  t = t ^ , the  shoreline  keeps  evolving  as  given  by 

equation  (11),  as  if  the  time  were  t - 0.38tj  . Then',  the  sediment  dis- 
charge, Q , bypassing  the  groin  is  equal  to  the  incoming  discharge 
minus  the  volume  of  sand  which  accumulates  per  unit  of  time. 


Q.t)  = Q„ 


tt  K ( t - D . 38t  j ) 


J 172  ; 


Q(t)  - Q0  / l - 


T 

1/2 

tana 

7iK(t  -0 . 38t  j ) 

0 

or  again 


Q(t ) = Q0  1 


0 . 638 


[(t/t1,  - 0 . 38  J 


In  dimensionless  terms,  the  following  values  are  obtained  for  equation 
f 1 <> ) (see  Fig.  5): 


' o 

0.189 

0.315 

0.397 

0.498 

0 . 605 


0 . 665 
0.703 


It  takes  a long  time  before  the  value  of  Q approaches  initial  dis- 
charge, , downdrift  of  the  groin. 
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r 


The  shoreline  may 

formation  about 

the  i 

given  time,  t , 

, and 

AD 

ft  - 0. 38t  j 

r 

be  deduced  at  any  time,  t , by  a homothctic  trans- 
oy  axis  from  the  knowledge  of  the  shoreline  at  a 
also  by  applying  the  simple  relationship  (sec  lig.3) 


AC! 


- 0.38t." 

1/2 

Z 1 

(17) 


The  theory  of  Pelnard-Considere  has  been  verified  in  laboratory  ex- 
periments with  fairly  good  accuracy.  The  steady-state  littoral  drift, 

Q , was  obtained  experimentally  from  preliminary  calibration  over  a 
straight  shoreline.  The  results  of  these  experiments  are  shown  in  Figs. 
6 and  7.  However,  the  shoreline  predicted  by  theory  is  not  expected 
to  be  valid  downdrift  of  the  groin  because  of  the  influence  of  wave 
diffraction  around  the  groin  tip.  Some  sand  begins  to  bypass  the  groin 
by  suspension  before  t = t^  (see  Fig.  5).  Also,  different  boundary 
conditions  apply  to  different  contour  lines  since  the  deeper  contour 
lines  reach  the  end  of  the  groin  before  the  contour  lines  which  arc  near 
the  shoreline,  which  implies  the  one-dimensional  theory  is  no  longer 
entirely  satisfactory. 

Subsequently,  Lepetit  (1972)  also  conducted  laboratory  experiments 
which  verify  the  results  of  a numerical  scheme  based  on  the  theory  of 

Pelnard-Considere.  He  used  the  law,  Q“  sin  a Vcosa"  . l.epctit's  ex- 
periments were  carried  out  with  a very  small  angle  between  wave  crest 
and  shoreline. 

1 . Refinement  and  extensions  of  the  Pelnard-Considere  Model. 

After  Pelnard-Considere ' s contribution,  the  mathematical  formula- 
tion of  shoreline  evolution  has  proceeded  at  a slow  pace.  The  first 
refinements  came  in  improving  the  longshore  transport  rate  (littoral 
drift)  formula,  in  particular,  modifying  the  expression  relating  sedi- 
ment transport  to  incident  wave  angle. 

Based  on  results  from  laboratory  experiments  performed  by  Sauvagc 

7a 

and  Vincent  (1954),  Larras  (1957)  introduced  the  function  f(a)  = sin  , 

also  used  by  Le  Mehaute  and  Brebncr  (1901).  New  theoretical  forms  of 
shoreline  evolution  are  determined  as  solutions  of  the  diffusion  equa- 

7ot 

tion.  Introduction  of  the  relationship  f(a)  = sin  -r—  instead  of  tana, 

4 

allows  obtention  of  solutions  valid  for  larger  wave  angles. 

Of  particular  interest  arc  the  cases  of  shoreline  undulations,  since 
assuming  linear  superposition,  any  form  of  shoreline  may  be  approximated 
by  a Fourier  scries.  The  solution  of  the  diffusion  equation  is  then  of 
the  form: 
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y 


cos 


Be 


-bt 


(X  - xo) 


( 18) 


b’ 


which  indicates  that  shoreline  undulations  tend  to  decay  exponentially 
and  disappear  with  time.  B defines  the  beach  undulation  amplitude  at 
time,  t = o , and  A is  related  to  the  wavelength,  I.  , of  this  undula- 
tion through  the  relationship: 


A = 


( 19) 


Shoreline  evolution  due  to  the  sudden  dumping  of  material  at  a given 
point  may  be  represented  by: 


>' 


7 


-x“/4Kt 


(20) 


Equation  (20)  gives  the  spreading  of  the  sand  along  the  shoreline  since 


the  integration 


/ 


ydx,  which  expresses  the  conservation  of  sedi 


ment  in  the  system,  is  a constant  (see  Fig.  8).  This  solution  was  also 
mentioned  by  Pelnard-Considcre. 


It  is  interesting  that  much  later,  Noda  (personal  communication, 
1974)  investigated  the  same  problem  by  taking  an  initial  condition  for 
sand  dumping. 

IY  = constant  when  |x| <X 
° when  |x| >X 


as  shown  on  Fig.  9.  Using  the  functional  relationship  now  commonly 
accepted,  f(a)  = sin  2a  , Noda  found  that  the  solution  to  the  diffusion 
equation  to  be: 
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Figure  8.  Spreading  of  sand  along  a shoreline  due  to  instantaneous 
dumping  at  a point. 


-X  X 


Figure  9.  Sand  dumping  along  a finite  stretch  of  beach 
(initial  condition). 
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liven  though  the  initial  condition  is  different  from  the  previous  one,  the 
solutions  tend  to  be  similar  as  time  increases  and  are,  therefore,  both 
applicable  to  the  problem  of  shoreline  sand  dumping. 

Also  of  interest  is  the  solution,  proposed  by  Larras  (1957),  of  a 
beach  equilibrium  shape  between  two  headlands  or  groins  described  by  the 
equation: 


3Q  - 

3s 


o , 


where  s is  the  distance  along  the  shoreline.  This'  indicates  no  sand 
transport  along  shoreline  configuration  and,  therefore,  yields  an 
equilibrium  to  obtain: 


7a 

ds  = L cos  — j da  (where  L is  a proportionality  constant), 
which  gives 


x = R 


T . 1 la  11  3a  1 

[Sln  — + “3  Sin  T J 


d f 1 lot  11  3a  I 

-R  [cos  — ~1  C0S  tJ 


(22) 


liquation  (22)  defines  a hypocycloidal  form  as  might  be  found  between  two 
headlands  (see  Fig.  10).  R is  a parameter  which  is  related  to  the 
relative  curvature  of  the  shoreline.  When  R ■*  <"  , a straight  shoreline 
solution  is  obtained. 

Another  family  of  solutions  was  given  by  Grijm  (1960,  1964).  In 
these  two  publications,  Grijm  used  the  most  commonly  accepted  expression 
for  dependence  of  longshore  transport  on  angle,  f(a)  = sin  2a  , and 
applied  the  theory  to  cases  where  the  angle  of  incidence,  a , is  not 
necessarily  small.  Subsequently,  he  established  the  kind  of  shoreline 
which  can  exist  mathematically  under  steady-state  conditions. 

Even  though  the  theoretical  approach  obeys  the  same  physical  assump- 
tion as  the  previous  theory  (except  for  the  allowable  range  for  the 
angle  of  incidence),  his  mathematical  formulation  is  not  as  simple.  The 
shoreline  is  defined  with  respect  to  a polar  coordinate  axis.  The  con- 
tinuity equation  is  solved  in  parametric  form,  which  is  integrated 
either  by  computer  or  by  graphical  methods.  Details  of  Grijm's  compu- 
tations are  not  available. 
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HYPOCYCLOID 


HEADLAND 


SAND 


Figure  10.  Equilibrium  profile  between  two  headlands. 
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The  main  interest  of  the  report  lies  in  the  results.  When  the  long- 
shore transport  rate  reaches  a maximum  value  (a  = 45°)  , the  shoreline 
tends  to  form  a cusp;  i.e.,  a cape  as  shown  in  Fig.  1). 

Also  of  interest  is  Grijm's  (1964)  mathematical  formulation  for 
different  forms  of  river  deltas  for  which  he  finds  two  possible  solu- 
tions, one  with  an  angle  of  wave  incidence  everywhere  less  than  45°,  and 
another  with  the  angle  of  incidence  greater  than  45°.  The  shoreline 
curvature  also  depends  upon  the  angle  a as  shown  on  Fig.  11.  The 
problem  remains  indefinite  since  it  is  unknown  which  solution  is  valid. 

The  formulation  of  Grijm  does  not  lend  conveniently  to  numerical 
adapt at  ion . 


Bakker* and  Edelman  (1964)  also  studied  the  form  of  river  deltas,  but 
instead  of  using  f(a)  = sin  2a  , as  Grijm,  they  used  the  linear  approxi- 
mation as  given  by  Pelnard-Considere;  i.e.,  f(a)  = k^  tana  for 
o < tana  < 1.23  . They  also  investigated  the  case  of  large  angle  of 
approach  using  the  function- 

k 7 

f(a)  = for  1.23  < tana  < °° 

tana 


Bakker  and  Edelman's  (1964)  solutions  are  similar  to  that  of  Grijm; 
however,  they  also  found  a periodic  solution  as  I.arras  (19S7)  did: 


= exp 

[ -dQ 

I da 

P 

II 

P 

o 

1 

2 

1 K~t 

D 

L 

J 

cos  kx 


(23) 


Equation  (23)  represents  a sinusoidal  shoreline  for  which  the  ampli- 
tude of  the  undulations  decreases  with  time  if  ^ is  positive  (i.e., 

da  ,q 

for  small  angles  of  wave  incidence),  but  increases  when  is  negative 

uct 

(i.e.,  for  large  angle  of  wave  incidence).  The  shoreline  is  thus  un- 
stable and  the  amplitude  of  the  undulations  increases.  It  can  be 
deduced  that  Grijm's  solution  for  large  angles  of  incidence  is  not 
natural ly  found,  since  they  are  unstable  and  will  be  destroyed  as  small 
perturbations  trigger  large  deviations. 


Bakker  (1968a)  implies  that  Grijm  did  not  discover  this  instability 
because  he  confined  himself  to  solutions  growing  linearily  with  t in 
all  directions,  while  the  exponential  solution  in  t also  exists. 

Komar  (1973)  also  applies  a numerical  scheme  based  on  the  Pelnard- 
Considere  approximation  to  the  problem  of  delta  growth.  He  found  shore- 
line shapes  identical  to  Grijm  in  the  case  of  a small  angle  of  approach. 

Irom  these  investigations,  it  is  remembered  that  the  Pelnard- 
Considere  approach  is  very  powerful  to  predict  shoreline  evolution  under 
small  angle  of  incidence.  But  under  large  angle  of  incidence,  instabili- 
ty of  the  shoreline  makes  it  very  difficult.  Furthermore,  the 
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(unstable) 

Figure  11.  Two  theoretical  forms  of  shoreline  equilibrium  of 
river  deltas. 
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f 


phenomenology  of  interaction  between  wave  and  shoreline  is  not  accurately 
defined  mathematically. 

2.  Example  of  Shoreline  Evolution. 

Because  of  its  importance,  an  example  application  of  the  theory  of 
shoreline  evolution  is  presented.  However,  the  example  is  slightly 
modified  to  account  for  the  generally  accepted  longshore  transport  rate 
formula: 

« ■ TE  »8  "b  cg  sin  2“b  - (2,) 

where 

Q = longshore  transport  rate  cubic  feet  per  second 
= wave  breaking  angle 
11^  = breaking  wave  height 
C = wave  group  velocity  at  breaking 

k = a constant  - 6.42  x 10  3 
pg  = specific  weight  of  seawater. 

For  the  case  of  a groin  perpendicular  to  shore,  consider  the  average 
beach  conditions: 

11,  = 5 feet 

b 

d,  = 6.4  feet 

0 

a = 5 

o 

D = 20  feet 

C = \/ gd,  = 14.4  feet  per  second 

g b 


^ 16  08  Hb  Cg  COS  2ab 

(6.42  x 10"3)  (64)  5 2 (14.4) 

8 x 20 


Substituting  into  equation  (10),  yields: 
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t 


1 


2 

H 71 

~> 

4k  tan“a 

o 


1.3  x 10 


-3 


2 

i days . 


In  tabular  form  for  various  groin  lengths. 


l ft 


50  100  200  500 


t days  3 


13 


52  325 


Check : 


For  i = 50  feet 


Area  Ox^y  = 1.56  x Area  oxy 


1.56 


_ 0.78  SL~  = 22,400  square  feet 

2 tana  tana 


Volume  = (Area  ox^y)  (D)  = 4.5  x 10  cubic  feet 
Q = — ^ tan  2a.  = 1.6  cubic  feet  per  second 


4.5  x 10  „ „ ,„5 

— — — pr = 2.8  x 10  seconds  = .•>  days. 


III.  THE  TWO- LINE  THEORY  OF  BARKER 


One  limitation  of  the  solutions  of  Pelnard-Considere  is  the  assump- 
tion of  parallel  depth  contours.  Bakker  (1968a)  realized  that  the  one- 
line  theory  of  Pelnard-Considere  and  its  subsequent  development  may,  at 
times,  lead  to  some  inaccuracy,  since  beach  slope  variations  along  the 
shore  were  not  considered.  Beach  slope  variations  with  respect  to  time 
(summer-winter  profiles)  arc  not  important  in  the  long-term  shoreline 
evolution.  Nevertheless,  if  an  adequate  onshore-offshore  profile 
response  model  was  available,  a suitable  mathematical  representation  of 
it  could  be  developed  (Dean,  1973;  Swart,  1974). 

Near  coastal  structures,  the  deviations  of  the  model  from  prototype 
conditions  can  be  considerable.  Pelnard-Considere  finds  that  the  accre- 
tion and  erosion  patterns  are  symmetrical  with  respect  to  the  groin  as 
shown  on  Fig.  12.  However,  in  reality,  the  updrift  profile  becomes 
steeper  than  the  equilibrium  profile  and  the  sand  moves  seaward.  The 
downdrift  profile  is  flatter  than  the  equilibrium  profile  and  the  sand 
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More  Reasonable  Approximation 


Figure  12.  Differences  on  shoreline  configuration  due  to 
onshore-offshore  transport  near  a groin  (from 
Bakker,  1968b). 


I 


r 


is  pushed  shoreward  by  the  waves.  To  reproduce  the  onshore-offshore 
movement  in  a mathematical  model,  it  is  necessary  to  schematize  the 
coast  by  two  or  more  contour  lines  instead  of  one. 

Bakker's  (1968b)  two-contour- 1 ine  model  is  not  easily  applied  to 
practical  engineering  problems  encountered  by  designers,  due  to  lack  of 
knowledge  about  onshore-offshore  transport.  However,  his  contribution 
toward  establishing  a realistic  mathematical  model  of  shoreline  evolu- 
tion is  of  sufficient  importance  to  deserve  detailed  review. 

Bakker  (1968b)  assumes  that  the  profile  is  divided  into  two  parts 
(Fig.  13).  The  upper  parts  extending  to  a depth,  Dj  , arc  affected  by 

the  groin,  the  part  below  1)^  extends  offshore  to  a depth  of  1) ^ + l)7 

which  is  the  assumed  practical  seaward  limit  of  material  movement. 

The  "equilibrium  distance",  w , is  defined  by  a distance  (y'  - y ) 

corresponding  to  an  equilibrium  profile  under  normal  conditions;  i.e., 
far  away  from  the  groins. 

The  onshore-offshore  transport  is  defined  by: 

Qy  = qy  Cy j - y;  + w)  . (25) 

where  q is  a proportionality  constant  (dimension  LT  ^).  When 

(y:  - y'  + w)  is  positive,  the  transport  is  seaward;  when  negative,  it 

is  shoreward.  q has  been  found  by  Bakker  for  a part  of  the  Dutch 
coast  equal  to  T to  10  meters  per  year  for  a depth  = 3 meters. 

Letting  y2  = y'  - w , then,  = qy  (y.  - y^  _ 

Now,  following  Pelnard-Considere;  i.e.,  developing  the  expression 
for  the  longshore  transport  rate  Q in  a Taylor  series  in  terms  of  a , 

Q = Q (a  -a)  +...  (26) 

o 9 a o 

which  gives  in  linear  approximations: 
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Defining 


(28) 


then. 


Q = 


q 


3y 

3x 


This  equation  is  now  applied  to  both  lines,  y ^ ( x ) and  y.,(x) 

3Xl 

Q1  " Qol  ' ql  3x 
3y2 

Q2  " Qo2  " q2  3x 
The  equation  of  continuity. 


(29) 


(30) 


=o 


(31) 


is  modified  by  the  term  due  to  onshore-offshore  transport  so  that 


- Qy 

9yi 

= D,  — — 

3x 

1 3t 

3Q2 

*Qy 

3y2 

= D.  — - 

3x 

2 3t 

(32) 


(33) 


Substituting  equations  (1),  (2),  and  (3)  for  Qj , Q-,,  Qy  gives: 

ay, 


2 

3 y 


11  3x2 


- - qy  (yi  ■ y2]  Di  at 


(34) 
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q>,  (y2  - yj) 


q2  — t 

" 3x" 


Adding  both  equations  yields: 


a i_z  + _U 

n 2 7 

3x  I) 


32(X1  " y2) 


in  which 


d = nTTl^  and  y = S’  (Diyi  + D2y2) 


q,  q2 

lor  simplicity,  Bakker  (1968b)  assumes  tt—  = j-—  , which  implies  that 

1 U2 

derivatives  of  the  littoral  drift  transport  with  respect  to 

do 

a:  are  proportional  to  depth  D . Then,  dividing 

a = a 

o 

equation  (6)  by  and  D7  respectively,  and  subtracting,  yield  : 


a2(yj 


3(yJ 


where  y*  = y.  - y_  , which  is  the  equation  for  the  offshore-onchore 
transport  q^y*  . It  is  interesting  that  the  offshore-onshore  transport 
is  independent  of  the  longshore  transport. 

Using  the  auxiliary  variable, 


ye  = y*  exp 


f 


f 


the  diffusion  equation  is  still  obtained: 


-> 


(40) 


Bakker  has  applied  his  theory  to  a number  of  idealized  cases, 
eluding  the  behavior  of  a sand  beach  near  a groin,  assuming 


in- 

141) 


The  boundary  conditions  are: 


a. 


b. 


Initial  condition  (t  = o)  : y^  = y?  = o for  o<x«=°  and  t = o , 
Then,  when  t>o: 

(1)  y,  = y,  = o for  x-*»  and  o<t<“>  (which  implies  an  equilibri- 

um profile) 

(2)  y = o for  x = o 


3y,  . xol  c 

(a)  1 1 = tana  for  x = o 

9x~  ql 


The  results  are  expressed  in  terms  of  lengthy  power  series,  and  are 
represented  graphically  in  Fig.  14. 


The  case  of  equilibrium  beach  profiles  between  groins  was  also 
investigated  by  Bakker  (1970). 

Despite  the  complex  refinement  of  the  two-line  theory,  as  initially 
developed  by  Bakker,  a number  of  phenomena  that  have  significant  in- 
fluence on  the  beach  profile  are  still  neglected.  Among  these  are: 

a.  The  influence  of  rip  current  near  the  groins  is  twofold:  rip 
currents  transport  material  from  beach  to  the  offshore  and  cause  wave 
refraction . 
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Figure  14.  Evolution  of  shoreline  and  offshore 
beach  limit  near  a groin 
(from  Bakker,  1968b). 


b.  The  influence  of  diffraction  on  the  leeward  side  of  groins  which 
has  an  effect  in  the  immediate  vicinity  of  the  shoreline. 

c.  The  effect  of  changing  wave  direction  caused  by  refraction  changes 
the  magnitude  of  longshore  transport  rate  and  the  boundary  conditions. 

d.  The  nonlinearity  in  th(p  transport  equation  is  of  minor  importance 
for  small  angles  of  incidence  (for  a>  45°  , the  coastline  becomes 
unstable  as  previously  mentioned). 

The  two-line  theory  has  been  verified  experimentally  (llulsbergen , 
Van  Bochove,  and  Bakker,  1976),  and  shows  a trend  identical  to  the  ex- 
perimental results.  There  are  some  differences  at  a small  scale  due  to 
secondary  currents,  breaking  wave  type,  changes  of  wave  height  due  to 
small  changes  in  morphology,  etc.  These,  however ,-%>re  short-term 
rather  than  long-term  evolution  phenomena. 


IV.  nil:  LI-'FliCT  OF  VVAVi:  DIFFRACTION 

The  effect  of  wave  diffraction  was  subsequently  taken  into  account 
by  Bakker  (1970).  Initially,  this  was  done  for  the  one-line  theory  of 
Pelnard-Considere  and  later  for  Bakker's  two-line  theory. 

Pelnard-Considere's  equat ions , 


Q • Qo,x.  - q(x)  g 


q(x) 


a 


a 


o 


(42) 


and 


9 y _ q(x)  3Q(x)  (43) 

3t  D 3x 


still  apply.  Q and  q are  now  functions  of  x , since  both  the 
incident  wave  height  and  angle  of  approach  vary  along  the  shore  with 
x , because  of  wave  diffraction. 

Inserting  the  expression  for  Q in  the  continuity  equation,  yields: 


9y 

at 


i_ 

i) 


3q  9y 
ax  ax 


, dQ 

1 o 
I)  dx 


144) 
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It  is  assumed  that  the  longshore  transport  rate,  Qq  , is  proportional 
to  the  angle  of  wave  incidence,  (u  - , and  the  Square  of  the  rela- 

X O X 

five  wave  height.  The  variation  of  wave  height  with  x is  given  by  the 
diffraction  theory  of  Putnam  and  Arthur  (1948).  The  modification  of 
wave  diffraction  by  wave  refraction  is  neglected. 


A similar  approach  has  been  proposed  by  Price,  Tomlinson,  and  Willis 


0 . 35 

(1972),  who  assume  that  Q = E sin  2a  , where  E is  the  trans- 

mitted energy  which  is  also  a function  of  x as  is  a (and  y is  the 

submerged  density  of  the  beach  material).  Price,  Tomlinson,  and  Willis 
then  obtain  the  one-line  theory  equation: 


0 . 35 


ys 


sin  2 a 


dE 

dx 


+ 2E  cos 


+ 


o 


(45) 


which  is  solved  numerically  with 


a 


a 


x 


-1 

tan 


3y 

d\ 


(40) 


Laboratory  experiments  were  performed  with  crushed  coal  by  Price, 
Tomlinson,  and  Willis  (1972).  The  theory  giving  the  effect  of  wave 
diffraction  was  verified  by  the  experiments  at  the  beginning  of  the 
test.  After  a 3-hour  test  which  may  correspond  to  a prototype  storm 
duration,  it  is  stated  that  the  wave  refraction  pattern  invalidates  the 
input  wave  data  and  a complex  boundary  condition  developed  at  the  up- 
drift  end  of  the  wave  basin. 

Bakker's  (1970)  consideration  of  wave  diffraction  has  been  included 
in  his  two-line  theory  where. 


*1 


q2 


3yl 
ql  3x 

• ^1 

(47) 

3y2 

II 

q2  3x~ 

(48) 

Neither  the  deepwater  line,  defined  by  y (x,t),  nor  q,  and  Q , is 

affected  by  diffraction.  Fig.  15  presents  typical  results  obtained  from 
this  theory  for  the  case  of  beach  evolution  near  a groin  and  between  two 
groins . 


* 
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Wave  incidence 


V.  SPIRAL  BLAU  II. S 


llooklike  beaches  (Fig.  1(>)  are  common  along  exposed  coasts  and  are 
formed  by  the  long-term  combined  effects  of  refraction  and  diffraction 
around  headlands.  Yasso  (1965)  discovered  that  the  planimetric  shape  of 
many  of  these  beaches  could  be  fitted  very  closely  by  a segment  of  log- 
arithmic spiral;  the  distance,  r , from  the  beach  to  the  center  of  the 
spiral  increasing  with  the  angle  0 according  to 


r = r 

o 


exp 


£o  cot 


■] 


(49) 


in  which  £ is  the  spiral  angle. 

Bremmer  (1970)  has  shown  the  logarithmic  spiral  to  give  an  excellent 
fit  for  the  profile  of  a recessed  beach  between  two  headlands. 

The  evolution  of  spiral  beaches  belongs  to  the  geographical  time- 
scale  domain  (Sylvester  and  llo,  1972).  However,  similar  evolution  has 
also  been  observed  over  smaller  time  scales  in  consonance  with  the 
definition  of  long-term  shoreline  evolution  adopted  in  this  study. 

So  far,  only  empirical  rather  than  theoretical  mathematical  repre- 
sentations of  spiral  beaches  are  available.  The  empirical  approach  has 
been  fruitful  in  providing  the  spiral  coefficients  £ as  function  of 
wave  angle,  a , with  the  headland  alinement  (Fig.  16)  (Sylvester  and 
Ho,  1972).  The  "indentation  ratio"  (depth  of  the  bay  to  width  of  open- 
ing) also  depends  upon  a and,  in  most  cases,  varies  between  0.3  and 
0.5  (Fig.  17). 

There  have  been  many  attempts  to  explain  this  peculiar  beach  forma- 
tion (Leblond,  1972;  Rea  and  Komar,  1975).  Leblond  assumed  that  the 
rate  of  sediment  transport  is  proportional  to  the  longshore  currents  as 
given  by  the  theory  of  Longuct-lliggi  ns  (1975).  He  also  assumed  that  the 
beach  profile  is  not  modified  by  erosion  or  accretion  so  that  the  con- 
tinuity equation  from  the  one- line  theory  can  be  used  in  a two-dimen- 
sional coordinate  system. 

Thus,  the  variation  in  longshore  current  intensity  with  wave  angle 
will  yield  the  rate  of  erosion  or  accretion. 

Difficulties  arise  in  expressing  this  variation  of  longshore  current 
in  areas  subjected  to  wave  diffraction.  Leblond  (1972)  points  out  that 
classical  wave  diffraction  theories  are  too  complicated  to  be  used  in 
his  theoretical  scheme.  Another  difficulty  arises  from  the  fact  that 
the  barrier  (headland)  is  not  thin  as  it  is  assumed  in  the  theory  of 
diffraction  of  Putnam  and  Arthur  (1948).  To  account  for  this  effect, 
Leblond  introduces  an  empirical  correction  coefficient  to  the  theory 
of  Putnam  and  Arthur  over  a two-dimensional  network.  The  results  of 


such  a complex  scheme,  which  is  plagued  with  numerical  instabilities, 
are  shown  in  Fig.  18.  Even  though  the  results  show  how  oblique  waves 
initiate  an  erosion  pattern  that  might  eventually  lead  to  the  formation 
of  hooklike  beaches,  they  do  not  show  that  the  beaches  represent  a good 
fit  to  segments  of  a logarithmic  spiral. 

Rea  and  Komar  (1975)  developed  an  approach  to  overcome  the  numeri- 
cal instability  encountered  by  Leblond.  They  combined  two  orthogonal, 
one-dimensional  arrays  as  shown  on  Fig.  19.  In  this  way,  deformation  of 
the  beach  can  proceed  in  two  directions  without  the  necessity  of  a two- 
dimensional  array.  The  wave  configuration  in  the  shadow  zone  was 
described  by  various  simple  empirical  functions  which  resulted  in  beach 
configurations  fairly  approximated  by  a logarithmic  spiral. 

- • 

The  main  interest  in  the  work  of  Rea  and  Komar  (1975)  is  that  they  show 
the  lack  of  sensitivity  of  the  shoreline  evolution  in  the  shadow  zone  to 
the  actual  pattern  of  incident  waves  used.  Also,  the  sensitivity  of  the 
beach  shape  to  the  energy  distribution  seems  to  be  small. 

VI.  PROTO  TYPE  APPLICATIONS 

The  application  of  mathematical  models  of  shoreline  evolution  to  pro- 
totype conditions  is  not  very  well  documented  in  the  literature.  It  is 
certain  that,  at  least  in  its  simplified  form  such  as  given  by  Pelnard- 
Considere,  the  method  has  been  used  by  practicing  engineers  and  designers. 

It  has  been  reported  in  unpublished  reports  but  very  little  has  appeared 
in  the  open  literature. 

Weggel  (1976)  has  formulated  a numerical  approach  to  coastal  process- 
es which  is  particularly  adapted  to  prototype  situations.  In  particular, 
it  includes: 

a.  A method  for  determining  the  water  depth  beyond  which  the  onshore- 
offshore  sediment  transport  is  negligible.  This  information  is  particu- 
larly useful  in  determining  the  quantity  D used  in  Pelnard-Considere' s 
theory  and  others.  It  is  also  useful  in  determining  the  effect  of  a 
change  of  sea  level.  Beach  profile  data  are  plotted  on  semilog  paper 
and  the  base  elevation  of  the  most  seaward  point  varied  until  an  approxi- 
mate straight  line  is  obtained  (see  Fig.  20).  He  found  D = 70  feet  at 
Pt.  Mugu,  California. 

b.  The  effect  of  a change  in  sea  level,  a situation  pertinent  to  the 
Great  Lakes,  is  also  taken  into  account  in  a way  proposed  by  Bruun 
(1962).  Using  the  principle  of  similarity  of  shoreline  profile,  the 
shoreline  recession  Ay  is  related  to  the  change  of  water  level  a by 
the  relationship  (Fig.  21): 


ELEVATION  ABOVE  SEAWARDMOST  POINT  (ft) 


Figure  20.  Semi logarithmic  profiles  (from  Weggel , 1976). 
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c.  A numerical  scheme  in  which  the  effect  of  wave  diffraction  could  be 
included. 

d.  A statistical  characterization  of  wave  climate  and  longshore  energy 
flux. 

Examples  of  recent  prototype  analysis  and  prediction  of  shoreline 
evolution  by  mathematical  modeling  are  Apalachicola  Bay  by  Miller  (1975) 
and  the  Oregon  coastline  by  Komar,  Lizarraga-Arciniega,  and  Terich  (1976) 
Both  studies  are  based  on  numerical  schemes  related  to  the  Pelnard- 
Considere  (one-line)  formulation. 


VII.  CONCLUSIONS 

There  are  two  methods  of  approach  to  the  problems  related  to  littoral 
processes.  The  first  one,  typified  by  the  previously  discussed  reports, 
consists  of  analyzing  global  effects.  The  method  essentially  based  on 
establishing  "coastal  constants"  for  a model  by  correlation  between 
long-term  evolution  and  wave  statistics  and  subsequently,  to  use  the 
model  for  predicting  future  effects.  It  appears  that  this  method  is  the 
most  promising  for  engineering  purposes  and  could  be  termed  the  macro- 
scopic view.  The  main  results  are  summarized  in  Table  2. 

The  second  approach,  the  microscopic  view  of  the  problem,  consists 
in  analyzing  sediment  transport,  step-by-step,  on  a rational  Newtonian 
approach,  starting  with  wave  motion,  threshold  velocity  for  sand  trans- 
port, equilibrium  profiles  of  beaches,  etc.,  until  the  individual  com- 
ponents can  be  combined  into  an  overall  model  to  predict  shoreline 
evolution.  The  second  method  or  scientific  approach  has  not  progressed 
to  the  point  where  it  can  be  applied  to  engineering  problems  in  the 
foreseeable  fut  ire. 

However,  much  progress  has  been  made  in  the  last  5 years  toward 
understanding  the  hydrodynamics  of  the  surf  zone  through  application  of 
the  "radiation-stress"  concept.  In  theory,  establishing  a reliable 
mathematical  model  of  surf  zone  circulation  should  permit  a determina- 
tion of  the  resulting  sediment  transport.  Practically,  however,  inter- 
action between  a movable  bed  and  the  surf  zone  circulation,  and  the 
inherent  instability  of  longshore  currents  limit  this  approach  to  the 
realm  of  research.  Among  the  problems  that  make  this  approach  difficult 
are  the  refraction  and  diffraction  of  water  waves,  uncertainty  in  pre- 
dicting rip  current  spacing,  and  the  effect  of  free  turbulence  generated 
by  breaking  waves  on  the  rate  of  sediment  suspension. 

Finally,  the  complexity  of  mathematical  formulation,  based  on  the 
radiat ion- stress  concept,  makes  it  difficult  to  use  as  a predictive  tool 
when  dealing  with  forcing  functions  expressed  by  statistical  multi- 
directional sea  spectra.  This  method  is  promising  in  explaining  local 
effects  (e.g.,  near  groins),  rhythmic  topography,  beach  cusps,  and  short- 
term evolution  due  to  unidirectional  sea  states.  All  these  effects  are 
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Table  2.  Summary  of  mathematical  models  for  shoreline  evolution 


Date 

1956 

1957 

1960 

1961 

1964 

1964 

196* 


1970 

1972 

1972 

1972 

1972 

197  J 

1975 

1976 

1*776 


Author 

Sediaent  transport 
alongshore 

Validity 

Sediaent  transport 
onshore-offshore 

Theoretical  developaents. 

Eapenaenta  1 
yen  fication 

Application  to  ideal 

cases 

Pelnard - 

Considers 

tan  a a»o 

o da 

Very  saall  angle 

No 

Diffusion  equation 
closed- fora  solution 

Laboratory 
with  puaice 

Groins 

Larras 

7a  av 

sin  t—  a»a  -*• 

4 o ai 

Sail  1 angle  (<2S°) 

No 

Diffusion  equation, 
closed-fora  solution 

No 

Groins  - sudden  duap 
sinusoidol  undulation 
equilihnua  shape 
between  groins 

Grija 

sin  2a 

Sin  all  angles.  In 
case  of  large  angle 
an  inconsistency  in 
the  assiaption 

H 

No 

Nonlinear  differential 
equation 

NO 

Foras  of  deltas 

Le  Mehaute , 
Rrebner 

7a 

SI"  4 

Diffusion  equation 
closed-fora  solution 

No 

Groins-sudden  diap 
sinusoidol  undulation, 
equilibria  shape 
between  groins 

Grije 

sin  2a  iaptied 

Saall  and  large 
angle 

NO 

Cylindrical  systea  of 
coordi nates- nuaer teal 
or  graphical  aethod 

NO 

Foras  of  deltas 

Bakker . 
Edelaan 

k tan  a 

* .i 

*2  tan  a 

f,  ’ *2  • '•« 

Saall  angle 
o<tana<l  .23 
large  angle 

1 . 23<tano • 

No 

Nonlinear  differential 
equations, 

closed-fora  solutions 

NO 

Foras  of  deltas 

Bakker 

tan  a 

Very  saall  angle 

Yes 

(two- line  theory) 

Systea  of  linear  differ- 
ential equations 

1)  power  series  solution 

2)  closed-fora  solution 

3)  closed-fora 

No 

Groins  and  coabma- 
tions  of  groins 

11  staple  groins 

2)  stationary  shore- 
lines 

3)  sand -wave  pro- 
pagation 

Bakker  . 
Breteler . 
Boos . 

sin  2a 

Saall  angle 

Nuaertca)  aethod 

No 

Groins 

Price. 

Toal inson , 
■litis 

sin  2a 

• ■ %-*-■% 

Saall  angle 

No 

Nuaerical  aethod  based 
on  Pelnard-Considere 
(1956) 

Laboratory 
with  crushed 
coal 

Groin 

Lepetit 

sina  cosa 

Saall  angle 

No 

Niaerlcal  aethod  based 
on  Pelnard-Considere 
(1956) 

Laboratory 
with  lakelite 

Groin  (updrift 
and  downdrift) 

Sylvester. 

Hi 

Eapirical  fit 

res 

Crenu 1 at ed - shaped 
bay  or  spiral 
beaches 

Leb 1 ond 

Proportional  to 
longshore  current 
f radiation-  stress) 

Saall  angle 

No 

Niaerlcal  aethod 

Spiral  beaches 

toner 

sin  2a 

Saall  angle 

No 

Niaerlcal  aethod 

No 

Growth  of  deltas, 
reorientation  of 
beaches  between 
two  headlands 

koaar 

sin  2a 

Any  angle 

Ntaertcal  aethod  based 
on  eapirical  aodel  re- 
fraction-diffraction 

Hooked  beaches 
(spiral  beaches) 

‘fcjlsbergen. 
vao  Bor hove . 
Bakker 

sin  2« 

Saall  angle 

Yes 

Application  of  the  two- 
line  theory  of  Bekker 

Laboratory 
with  dune  sand 

Groin 

■**gel 

sin  2a 

Saall  angle 

Yes 

Natheaatical  and 
nuaerical  foraulation 

No 

Groin 
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Table  2.  Summary  of  mathematical  models  for  shoreline  evolution . --continued 


Variation  of  beach 
slope  taken  into 
account 

M«>d  1 1 ic.it  ion 
Effect  of  by 

Diffraction  refraction 

Variable 

wave  Variation  of 

direction  sea  level 

Effect  of 
rip  current s 

Application  to 
prototype  cases 

Main  Conclusions 

Mo 

(one-line  theory) 

NO 

No 

No 

No 

No 

No 

First  significant  milestone  of  intro- 
duction of  mathematical  model!' g to 
the  study  of  shoreline  evolution 

Mo 

NO 

No 

No 

No 

No 

No 

Eitens ion  of  the  method  of  Pelnard 
Consider*  to  other  idealized  cases 

Mo 

NO 

No 

No 

No 

No 

No 

Two  forms  of  solution 

One  with  concave  shoreline  (small  angle) 
One  with  conves  shoreline  (large  angle  1 

No 

NO 

No 

No 

No 

No 

No 

Same  as  Lerrss  (I9S7) 

No 

NO 

No 

No 

No 

No 

No 

Two  forms  of  solutions  as  in  I960 
applied  to  a mmher  of  deltas, 
idealized  cases 

Mo 

NO 

No 

No 

No 

No 

No 

Instability  of  shoreline  under  large 
angle  (even  a straight  shoreline) 

Implicitly  (through 
the  two- line  theory) 

NO 

No 

No 

No 

No 

No 

The  most  significant  contribution  since 
I9S6  demonstrating  the  influence  of 
beach  slope 

Implicitly 

Ye* 

(periodic  wave) 
(constant  depth) 

No 

Yes 

No 

No 

No 

Influence  of  diffraction  changing  wave 
condition  leading  to  stable  shoreline 
near  groin 

Mo 

Yes,  but  not  cam- 
pletely  formulated, 
not  applied 

No 

No 

No 

No 

No 

Espenmental  verification  for  more 
diffracting  waves,  i.e.,  updnft, 
fairly  satisfactory 

No 

No 

No 

No 

No 

No 

No 

food  experimental  verification 

Combined  effect, 

refraction- 

diffraction 

Yes 

No 

No 

Fit  with 
prototype 
cases 

Combined  effect  of  refraction  i dif- 
fraction in  protected  areas,  affected 
also  to  geological  time  evolution 

Mo 

Yes,  but  unsuccess- 
fully / manor  i cal 
instability) 

Yes 

No 

No 

No 

No 

Qualitative  fit  only,  unsuccessful,  two 
large  distances  between  oats  point . 
Complex! t\  of  combined  refraction- 
diffraction  effect . 

No 

NO 

No 

No 

No 

No 

No 

Nimwrical  application  of  Pelnard- 
Considere 

Combined  re frac- 
tion and  diffrac- 
tion 

Empirical  development 

Implicitly 

No 

No 

No 

No 

No 

No 

Fairly  good  experimental  verification 
except  near  groin 

Yes 

(in  principle) 

Yes 

(in  principle) 

No 

Yes 

(statistically) 

Yes 

No 

No 

Aimed  at  investigating  real  cases 
on  Greet  Lakes 
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superimposed  on  the  long-term  evolution  for  which  an  analysis  can  be 
done  independently. 
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Among  the  significant  recent  reports  leading  toward  understanding  of 
surf  zone  circulation  and  related  bottom  topography  are:  Bowen  and 
Inman  (1969)  who  advocate  the  presence  of  edge  waves  as  a cause  of  rip 
currents  and  beach  cusps;  llino  (1974)  who  states  that  rip  currents  are 
the  result  of  mobility  of  the  sand  bed  and  hydrodynamic  instability; 

Sonu  (1972)  and  Noda  (1972)  demonstrated  that  a perturbation  on  bottom 
topography  causing  waves  to  refract  and  have  varying  intensity  along  the 
shore  induces  a variation  in  radiation  stress  which. in  turn  enhances  rip 
currents;  finally,  Liu  and  Mei  (1976)  applied  the  radiation-stress 
concept  to  a groin  perpendicular  to  shore  and  to  an  offshore  breakwater. 

These  investigations  offer  at  least  partial  answers  to  a number  of 
important  problems,  important  in  understanding  shoreline  processes.  It 
definitely  indicates  that  the  radiation-stress  approach  holds  the  poten- 
tial key  to  understanding  many  types  of  nearshore  currents,  heretofore 
unexplored.  It  is  also  evident  that  the  study  of  surf  zone  hydrodynamics 
will  rapidly  reach  a plateau  if  sand-water  interaction  problems  are  not 
mastered,  and  at  this  stage,  these  can  only  be  considered  empirically. 
Determinism  leaves  off  with  the  inception  of  turbulence. 


Even  though  the  dynamics  of  nearshore  currents  hold  the  key  to 
understanding  of  beach  processes,  application  of  the  methodology  based 
on  radiation  stress  to  investigate  shoreline  evolution  mathematically  is 
still  beyond  the  state-of-the-art. 

Both  approaches  could  be  pursued  in  parallel  and  the  results  of  the 
scientific  approach  could  slowly  be  incorporated  into  a practical 
engineering  model. 

Conclusions  based  on  the  literature  survey,  as  summarized  in 
Table  2,  are: 

a.  There  is  sufficient  laboratory  verification  to  give  credibility  to 
a mathematical  approach  to  t lie  study  of  shoreline  evolution  for  small 
angles  of  wave  approach. 

b.  For  large  angles  of  incidence,  there  is  a lesser  chance  at  arriving 
at  a successful  formulation  since  shorelines  are  then  unstable  and  the 
resulting  shoreline  evolution  could  not  be  predicted  without  t1'  • 
initiation  of  more  basic  research  beyond  the  present  state  of  knowledge. 

c.  Even  though  no  field  measurements  subsecpicnt  to  mathematical  pre- 
dictions have  been  found  in  the  literature,  many  practicing  engineers 
have  applied  the  theory  of  Pelnard-Considcrc  (1956)  to  predict  shore 
evolution  by  taking  into  account  variable  wave  climate.  The  method  is 
easy  to  apply  and  provides  valuable  information. 
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<J.  Engineering  applications  to  prototype  cases  based  on  more  sophisti- 
cated approaches  such  as  given  by  the  two- line  theory  of  Bakker  (1968b) 
are  not  known.  These  more  sophisticated  approaches  can  be  currently 
considered  as  belonging  to  the  realm  of  research  rather  than  of  engineer- 
ing practice. 

e.  - Local  effects,  diffraction,  rip  currents,  wave  refraction  and  inter- 
action between  these  effects  are,  at  present,  still  not  so  conveniently 
formulated  to  be  used  by  practicing  engineers.  Introduction  of  these 
effects,  if  and  when  important  in  the  mathematical  formulation,  is 
feasible  but  will  require  further  investigation. 

f.  A simple  numerical  scheme  that  can  be  used  by  design  engineers  and 
planners  and  which  includes  theoretical  or  empirically  all  important 
effects  could  be  developed.  Effects  that  should  be  included  in  the 
mathematical  model  are  wave  diffraction,  loss  of  sand  by  rip  currents 
along  groins,  sea  (lake)  level  variation,  and  beach  slope  variation 
near  groins. 

g.  The  introduction  of  the  concept  of  radiation  stress  in  the  mathe- 
matical formulation  is  not  recommended  at  this  time,  but  research 
related  to  this  approach  should  be  pursued  in  view  of  the  eventual 
input  that  subsequent  results  could  have  on  then  existing  operational 
mathematical  models. 
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